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We present a hydro-dynamic theory of polar active smectics, for systems both with and without 
number conservation. For the latter, we find quasi long-ranged smectic order in d = 2 and long- 
ranged smectic order in d = 3. In d — 2 there is a Kosterlitz-Thouless type phase transition horn 
the smectic phase to the ordered fluid phase driven by increasing the noise strength. For the number 
conserving case, we find that giant number fluctuations are greatly suppressed by the smectic order; 
that smectic order is long-ranged in d — 3; and that nonlinear effects become important in d = 2. 

PACS numbers: 05.65.+b, 64.70.qj, 87.18.Gh 



Active matter [l[ can exhibit a richer variety of or- 
dered phases than for equilibrium systems, since symme- 
try differences which have little or no effect in equilibrium 
can have radical effects on active systems. For example, 
in equilibrium, systems with "polar" orientational order 
(e.g., ferromagnets) and those with apolar orientational 
order (i.e., nematics)@ have identical scaling of their ori- 
entational fluctuations. In contrast, the active polar or- 
der of a coherently moving flockH fluctuates far less than 
that in active apolar orientationally ordered systems (i.e., 
"active nematics" Q ) . 

In this paper we formulate the hydrodynamic theory of 
polar active smectic systems, by which we mean systems 
that spontaneously form uniformly spaced moving lay- 
ers, and find that they differ considerably from their ap- 
olar (i.e., non- moving) analogs Examples of such "ac- 
tive polar smectics" include propagating waves in chem- 
ical reaction-diffusion systems [y], and "flocks" @ of ac- 
tive particles forming uniformly spaced parallel liquid- 
like layers (density waves). Layers are ubiquitous [7| as 
the order-disorder transition is approached in, e.g., the 
"Vicsek model" @ of nocking, although they may not 
reach a uniform steady state spacing. 

We restrict ourselves here to "active smectics A" , 
meaning phases with the average particle velocity along 
the mean layer normal. "Active smectics C" @, with other 
relative orientations of particle velocity and the layers, 
will be considered elsewhere [To| . We consider two cases: 
first, with no conserved quantities, and second, with only 
particle number conserved. In neither case is momentum 
conserved. 

For the non-number conserving case (hereafter the 
"Malthusian" casell|), we find that the active polar 



smectic phase is stable over a finite range of parameters; 
that is, it is both hydrodynamically stable, and robust 
against noise. In particular, we find quasi-long-ranged 
smectic order in d = 2, and long-ranged order in d = 3. 
In contrast, thermal fluctuations completely destabilize 
the equilibrium smectic phase in d = 2 [12j , and allow 
only quasi-long-ranged equilibrium smectic order in d = 3 

m. 



More specifically, in d = 2 

(P(f,t)i/>(f',t)) cx \r-f'\-\ 



(1) 



where ij) is the complex smectic order parameter, defined, 
as in equilibrium smectics via : 



p(f,t) = p + i>(r,t)e iq °* 



(2) 



where p is the number density (or, more generally, the 
spatially modulated scalar field in the problem, such as 
chemical concentration in reaction-diffusion systems), po 
its mean, and qo = 2n/a, with a the distance between 
neighboring layers. Here we've defined the average plane 
of the layers as the _L plane, and the normal to this plane 
as the z axis. Equations (JTJ and (J2J imply quasi-sharp 
Bragg peaks in light or X-ray scattering (or, equivalently, 
the numerical Fourier transform of density correlations); 
that is: 



I n (q) oc (\p(q,t)\ 2 ) oc (q z - nq Q ) 2 + 79 



(3) 



where 7 is an O(l) constant, n is an integer denoting the 
order of the Bragg peak, and 77 is non- universal (i.e., it 
varies from system to system). 

For a finite system, this divergence is cut off for \Sq\ ~ 
l/L, where L is the spatial linear extent of the system, 
and Sq = q — nq^z. This implies the nth peak will have 
a finite height which scales with L like L 2 ~ n 17 . 

In d — 3, the long-ranged nature of the smectic order 
implies sharp (i.e., 8- function) Bragg peaks. In a finite 
system, the height of these peaks scales linearly with sys- 
tem volume L 3 . 

We predict that with increasing noise the active smec- 
tic phase undergoes a dynamical phase transition into 
the fluid, polar ordered phase treated in much prior 
workQ. This transition is in the equilibrium XY uni- 
versality class 3, which in d = 2 is of the the Kosterlitz- 
Thouless tvpe[14j|. The phase diagram in the parameter 
space of our model is illustrated in Fig. [2] 

In the case where the number of the particles is con- 
served, we find long-ranged smectic order in d — 3. In 
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d = 2, the linearized version of the full hydrodynamic 
theory predicts quasi-long-ranged smectic order; how- 
ever, there arc marginal non-lincaritics which may inval- 
idate this conclusion. We'll investigate this in a future 
publication [l5j. 

In neither case are there giant number fluctuations in 
d = 3; nor are there are any in d = 2 in the Malthu- 
sian case. The linearized hydrodynamic theory predicts 
none in d = 2 for the number conserving case either; but 
the aforementioned non-linear terms could change this as 
well. 

We'll now outline the derivation of these results, start- 
ing with the case without number conservation. Then the 
only important hydrodynamic variable is the displace- 
ment u(f,t) of the layers along z, which is proportional 
to the phase of the the smectic order parameter ijj{f,t): 

Symmetry considerations (specifically, translation and 
rotation invariance) require that it's equation of motion, 
to lowest order in a gradient expansion, take the form: 



d t u 



v Q - 2\ ± d z u + (vj_V 2 



+a 2 (d zU ) 2 + \x | v±u | 2 +/, 



(4) 



where / is a Gaussian, zero-mean, white noise with vari- 
ance {f(f,t)f(r ',*')) = 2A5 d (f-f')5(t-t'). Rotation 
invariance forces the coefficient of the d z u term to be ex- 
actly —2 times that of the | Vj_u | 2 term, because only 
the combination d z u — h | \7±u | 2 is unchanged by a 
uniform rotation of the smectic layersQ. 

The term with coefficient u± in (j4j) is forbidden by 
rotation-invariancc of the free energy in equilibrium. It is, 
however, permitted herefl7l[l8j simply because rotation- 
invariance at the level of the equation of motion, which 
is all one can demand in an active system, does not rule 
it out. Its physical content is that layer curvature pro- 
duces a local vectorial asymmetry which must modify the 
directed motion of the layers as this is a driven system. 

In contrast to apolar active smecticsJH, Eq. @ is 
not invariant under the simultaneous transformation u — > 
—u, z — > —z, due to the lack (by definition) of up-down 
symmetry in polar smectics. 

To simplify Eq. (gj), we introduce another field variable 

u 1 = u — VQt (5) 

and another coordinate system t' =t,z' = z — 2A±t, r_/ = 
f±. In terms of these Eq. ((4| becomes 

dt-u = v±Vl,u' + v z dt,u' + \±\V±>u'\ 2 + \ z (d z >u'f 



+/, 



(6) 



In Eq. ([5]), v_\_. z must be positive for the smectic state to 
be dynamically stable. Either sign of Aj_ jZ can be stable; 
indeed, their signs need not be the same. 

Eq. (JSJ has exactly the same form as the anisotropic 
KPZ equation [l9j]. However, there is a crucial difference. 
The original KPZ equation[2(| describes the hydrody- 
namics of crystal growth, and the hydrodynamic variable 



is h, the height of a d dimensional surface. Clearly, states 
with different heights h are always physically distinguish- 
able. However, for smectics, the state is periodic in u' 
with period a, the spacing between neighboring smec- 
tic layers. This allows for the existence of topologically 
stable dislocations, which can unbind, thereby "melting" 
(i.e., disordering) the smectic, in analogy to such "dislo- 
cation mediated melting" in a variety of translationally 
ordered equilibrium systems fl2l [2lj . In d = 2, this is 
the aforementioned Kosterlitz-Thouless phase transition, 
which is absent in the anisotropic KPZ equation. 

Simple power counting shows that the nonlinear terms 
in Eq. © are irrelevant in d = 3; hence, the linear theory 
is valid. A straightforward calculation then shows that 
(\u'(q, t) | 2 ) cx 1 /q 2 for all directions of wavevector q. This 
in turn implies that the real space fluctuation (\u'(r, t)\ 2 } 
is finite as system size L — > oo, which implies long-ranged 
smectic order (22^. 

In d = 2 the nonlinear terms in Eq. ([5]) become 
marginal, and a dynamical renormalization group (RG) 
analysis is needed. This has already been done for the 
d = 2 crystal growth problem the resulting RG re- 
cursion relations are: 



dvx 
d\± z 



- 2 + 5 fc« 1 - r > 



= (x + z-2)X^ z 



(7) 
(8) 



dA 

U 

d{v z /v\_) 



-2 X + z - 2 + -f - (3r 2 + 2r + 3) A, (9) 
647T v '\ 



= -JLi^ (l _r 2 ) 

di 32tt v± v ' 



(10) 



where \ and z are the rescaling exponents of u' and t, 



(i.e., u' -> u'e xt , t' -> t'e zl ), V = = ^\ 



we have chosen to rescale lengths isotropically. Eqns. 
([71 fTU|) imply a set of closed flow equations in T-g space: 



HI 
dg 



IJL 

32tt 
32tt 
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(r 2 



■4T- 1 



(11) 
(12) 



The ratio of these two equations yields a separable 
ODE for r as a function of g. Inserting the resultant 
solution for r(£/) into (TT2"|) yields an integrable equation 
for g(l); inserting that g(l) into pT|) yields a solvable 
equation for T(£). We thereby find 



(i + iy 



ro.9o 

i-r 
i -r 



i-r 



i + r 

2 



r n 



(i + r r 



l + r 
i + r 



(13) 



(14) 



where go and T denote respectively the "bare" values of 
g and T (i.e., their values at £ = 0). This solution implies 
a stable fixed point for Tq < 0: as I — > oo, T(£) —> — 1, 
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o 41-2 



FIG. 1: The RG flow in the T-g parameter space for active 
smectics in d = 2. For T < and g > 0, all flow lines go to a 
stable fixed point ( — 1,0); for T > and g > 0, all flow lines 
go to infinity. 



1 + T(£) oc l/y/t, and g{£) oc l/l The RG flow loci inT-g 
space are illustrated in Fig. Q] 

We will now use the trajectory integral matching 
method Q to compute (\u'(q, t)\ 2 }, which determines the 
presence or absence of smectic order |22|. We restricted 
this calculation to the case r < 0, since, as we'll see, 
only then is a stable smectic phase possible. Performing 
this standard procedure, we obtain 



v ± {£*) 



(15) 



where £* = InA/g, with A the ultraviolet cutoff. For 
small q (<C A), £* >• 1; in this limit, we find 



where l\ is some fixed value of the renormalization 
group "time" £ at which the large £ approximations be- 
come valid, and C is a constant that we could express 
in terms of Fo and go using Eqns. (fTB"]) and (|14p . if we 
cared. 

Note that the integral over £' in this expression con- 
verges as £ — > oo; hence, (£* — > oo) — > C'e~ 2xt , where 
C" is a finite, non-zero constant. This convergence makes 
the scaling of u' correlations the same as that predicted 
by the linear theory, as we'll now show. 

Using Eq. §T§§ in Eq. ([15]) gives 



(W(q,t)\ 2 ) = 



C" 



9 



hi 



(17) 



Now note that, since V 



— 1 as 



°°> 177(F) -> 



recursion relation 



ratio 



the last equality following because the 
does not renormalize, as can be seen from the 
for the A's. Thus we have, finally, 



(\u'(.q,t)\ 2 ) = 



a 



(18) 



— (£*) = exp 



-2 X (£*-h) + C 
A 



df 



£13/2 



(16) 



which clearly scales as l/q 2 for all directions of wavevec- 
tor q. This scaling implies "logarithmic roughness" of 
the smectic layers in d = 2: that is, 



u'(r,t) - u'{r',t) 



which implies Eq. (1) via [22| 



(ip* (r,t)ip(r',t)) oc exp 



1 2 
- 2 % 



A° 



u(f, t) — u(r' , t) 



hi 



(rj 



2" 









(r±-r' ± y 



AS 



(z-z'Y 



A3. 



A° 



(z-z'r 



(19) 



oc|r-r'|-", (20) 



with rj = C'qQ^/\ Aj_/A9 |/47r; i.e., quasi-long-ranged smectic order[14|]. This is the first real demonstration 
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that the roughness of the anistropic KPZ equation is only 
logarithmic in d = 2; earlier arguments did not address 
the possibility that the integral in equation (fTB")) could 
fail to converge as I — > oo, thereby changing the scaling 
from that predicted by the linearized theory (as happens 
in d = 3 equilibrium smecticsfHJ). 

When there is no number conservation, as in the model 
we are considering here, we will show later that, in both 
d = 2 and d = 3, the spatial Fourier transform of the 
number density p is given by 



{\p{q,t)\) = Ciql{\u{q,t)Y 



Co 



(21) 



where C\ and C2 are constants. Since this remains finite 
as q — » 0, there are no giant number fluctuations [25j . 

Now we discuss the stability of the smectic phase in 
d = 2. Doing an anisotropic rescaling r'j_ 



\J Ux/uzZ' ', expressing it' in terms of = 2ttu' /a, and 
ignoring the nonlinear terms, we can write Eq. as 



d t ,e = u A y' n e + f, 

where the statistics of /' are given by 

(f (r ", t')f (6,0)) = K v ± 5(r") 5(f) 
with k = A (27r/a) 2 j *Jv~\ 



(22) 



(23) 



_i/ z . Eq. (f22j) is identical to the 
simplest relaxational equation of motion for an equilib- 
rium XY model with 9 being the angle of the magneti- 
zation, and kv^ = 2ksT. This mapping implies a dislo- 
cation unbinding phase transition [2l| in active smcctics 
in d = 2, when k = 7r. For k > 7r, the system is in the or- 
dered fluid phase; for k < tt, the system is in the smectic 
phase. 

The effect of the nonlinear terms can be included sim- 
ply by replacing the bare n with its (finite) renormalized 
value. This implies that the transition occurs when 



lim k(£) 



(24) 



Note that Eq. (|2~4"|) is valid only if the period of the 
rescaled system on u' is kept fixed at a, which requires 
no v! rescaling (i.e., X = 0). Using our earlier solution to 
the recursion relations it is straightforward to show that 



lim k(£) -> -- — — k - 

1^00 w 4r 



Using this result in Eq. (jM)) we obtain the phase bound- 
ary in terms of kq and Tq for Tq < 



Kq 



_JttTo_ 

(i-r ) 5 



(26) 



This phase diagram in the Ko-Eo parameter space is il- 
lustrated in Fig. 

Now we turn to the case where the number of particles 
is conserved. In this case the fluctuation 8p = p — pa of 
the density p about its mean value po becomes another 











Active polar 


ordered fluid 




K 


Active smeqtic phase \ 





(25) ©. 



FIG. 2: Phase diagram in the kq-Fo parameter space. 

important hydrodynamical variable. Number conserva- 
tion implies dtSp = —V -j, where j is the number density 
current. Based on symmetry arguments, a gradient ex- 
pansion of j, keeping only "relevant" terms is given by 
j = jl + Jn l 1 where the linear piece 

Jl = - [jo + v p 5p + D z d z Sp + v z pu d z u 
+ ((c ± - w)V\ + c z d 2 z ) u]z 
-D ± V ± Sp - v^Vxn - wV ± d z it - f p , (27) 

where f p is a Gaussian noise with statistics 

(fpi(r, t)f pj (6, 0)) = (A 2 <% + Al^) S(r)S(t) , (28) 
while the non-linear piece is given by 

3nl = -[A^plV^ul 2 + X zp (d z uf + v p 5pd z u + g p 6p 2 ]z 
+g u d z uV ±u + v p 6pV ±u . (29) 

Similar symmetry arguments and gradient expansions 
give the equation of motion for u, again keeping all rele- 
vant terms, 

d t u = vq + v u d z u + v up 6p + v z o^u + v\p\u + v p d z Sp 
+A ± |V ±U | 2 + A Z (8 z uf + g8p 2 
+g c Spd z u + f u , (30) 

where the noise f u has the same statistics as / in Eq. 



If we neglect the non-linear terms in j and (1301) . 
a straightforward calculation shows that (\u(q,t)\ ) ~ 
1/g 2 , which implies quasi- long-ranged smectic order in 
d = 2 and long-ranged order in d = 3. We also find that 
(\Sp(q, t)\ 2 ) goes to a finite value as q —> 0, which implies 
no giant number fluctuations in either d = 2 or d = 3. 

Simple power counting shows that the nonlinear terms 
in j and (|30[) are irrelevant in d = 3, in the RG sense. 
Hence, these linear results should apply in d = 3, at 
least in systems with sufficiently small non-linearities [26j . 
In d = 2, similar power counting shows that all of the 
nonlinear terms in j and (|30[) become marginal, and, 
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hence, could potentially change the behavior at long 
wavelengths. 

Supplementing the continuity equation for Sp with a 
source term to reflect the tendency of birth and death 
to restore the local population density to its equilibrium 
value, and an additional, non-number-conserving noise 
reflecting statistical fluctuations in the local birth and 
death rate, as has been done[ll[ for flocks with polar ori- 
cntational order, enables us to analyze number fluctua- 
tions in the Malthusian well. Dropping irrelevant 
terms, we obtain the equation of motion: 



d t Sp = ad z u + v^V\u + v z pu 8 2 z u - S P /t + h-d , (31) 

where r is the characteristic relaxation time for density 
fluctuations to relax away due to birth and deathfllj], the 
a term reflects the fact that symmetry allows the local 
birth and death rate to depend on the local layer spac- 
ing, and fb-d is the aforementioned noise in the birth 
and death rate. We take fb-d to be zero mean Gaus- 
sian white noise, with statistics (fb-d(f, t)fb-d(f' ',t')) = 
2A b - d S d (r-r>)S(t~t>). 

Fourier transforming equation pip in space, and solv- 
ing the resultant linear stochastic ODE for the correla- 



tions of Sp gives 27], to leading order in q, equation (|21[) . 

■ . i 9 9 n ^ A 



with C\ = r 2 a 2 and Ci = tA 



b-d- 



In conclusion, we have developed the hydrodynamic 
theory of apolar active smectic in both d = 2 and d = 3. 
We considered both the case in which the number of the 
particles is conserved, and that in which it is not. We 
have made various experimental predictions which can be 
tested. For the number conserving case an RG analysis 
of the nonlinear terms remains to be done. 
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